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Abstract :Using a D meson cloud model we calculate the squared charm radius of the 
nucleon . The ratio between this squared radius and the ordinary baryon squared radius 
is identified with the probability of "seeing" the intrinsic charm component of the nucleon. 
Our estimate is compatible with those used to successfully describe the charm production 
phenomenology. 



In the early eighties there was a hope to understand charm production solely in terms of 
perturbative QCD. Inspite of all uncertainties in defining the scale , it would be in any case 
of the order of a few GeV and therefore the coupling constant would be smaller than one. 
As more and more data became available it became clear that perturbative QCD alone was 
not enough to properly account for the measured differential cross sections. Higher order 
corrections did not make the situation any better. The main problem was that the produced 
charmed particles were too fast. In other words, there was a remarkable excess of particles 
with large Feynman x { xp )■ 

Already at that time, the idea was advanced , that the hadron wave function contains a 
charm component even before undergoing a collision. This component is originated in higher 
twist QCD interactions inside the hadron. The so called "intrinsic" charmed pairs produced 
by these interactions have nothing to do with usual sea quark pairs. The crucial difference 
between them is that the intrinsic charm is part of the valence system and therefore very 
fast in contrast to the sea charm, which is slow. During the last years, an intrinsic charm 
component was added to the perturbative QCD component in a quantitative and systematic 
way PI . As a result , a very good description of data was achieved. In order to obtain 
such good agreement with experimental data the crucial point was the normalization of the 
intrinsic charm component , of the hadron + nucleon — c — c X cross section . The 
quantity cTjc is related to the probability of observing the intrinsic charm component of the 
hadron, Pj^ 0]- It is very difficult (if not impossible) to calculate this quantity from first 
principles. It was estimated from a phenomenological analysis to be less than 1% [Q. In 
fact. Pic = 0.3% seems to be the best value to describe recent data on charm production [0. 

A very important question is, of course, whether this 1% of intrinsic charm can be 
supported by any model based calculation. In ref. 0, such a calculation was done using 
the MIT bag model. It was found that the probability of finding a five-quark component 
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\uudcc) configuration bound within the nucleon bag is of 1 or 2%, in good agreement with 
the above mentioned phenomenological estimate. 

In this note we calculate Pic using an approach, which is completely different and in- 
dependent from that used in ref. [1-5] and can therefore be used as a cross-check to those 
estimates. 

The existence of intrinsic charm is here associated with low momentum components 
of a virtual c-c pair in the nucleon. At low momentum scales , the virtual pair lives a 
sufficiently long time to permit the formation of charm hadronic components of the nucleon 
wave function. The same argument is true for the strange matrix elements. 

Generally speaking, we can say that the proton is a fluctuating object, being sometimes 
a neutron plus a pion, sometimes a strange hyperon plus kaon and so on. It can be any 
combination of virtual hadrons possessing the right quantum numbers. In particular, if 
charmed pairs pre-exist inside the nucleon, it can oscilate into a charmed hyperon plus a D 
meson, as e.g., by the process 



We calculate the intrinsic charm contribution to the matrix element (A^|c7^c|A^) arising 
from this virtual D meson cloud. The idea that intrinsic quark contributions to nucleon 
matrix elements can be given by meson clouds is not new. It was used in refs. 0, 0) Hj HI to 
estimate the intrinsic strangeness content of the nucleon and it was suggested in ref [Q as a 
picture to understand the existence of intrinsic charm in the nucleon. 

As in ref. [Q, we compute the D meson loops using an effective meson-nucleon vertex 
characterized by a monopole form factor 



and we introduce "seagull" terms in order to satisfy the Ward-Takahashi (WT) identity. In 
eq. (2) m is the meson mass and A is the effective cut-off. The inclusion of the meson-nucleon 



p Ac + D ^ p . 




- A^ 
A;2 - A2 ' 



(2) 
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form-factors is important to properly take into account the underlying nucleon structure and 
its spatial extension. As shown in ref. P], when the sub-structure of the nucleon is considered, 
it is the size of the proton, rather than the masses involved in the loop, that determines the 
effective momentum cut-off. We expect therefore the effective cut-off in the D meson-nucleon 
form factor to be approximately the same used in the pion-nucleon or kaon-nucleon form 
factors. 

The pseudoscalar meson-baryon coupling for extended hadrons is schematically given by 

Cbbm = -tgBBA4^l5'^F{-d^)(l) , (3) 

where and (p are baryon and meson fields respectively, F{k'^) is the form factor at the 
meson-baryon vertices and k is the momentum of the meson. The fact that the nucleon-D- 
Ac coupling constant is not known is not important here because we are mostly interested 
in arriving at some upper limit to the intrinsic charm content of the nucleon and not at 
definitive numerical predictions. Accordingly we will use the pion-nucleon coupling constant 
as an upper hmit to the nucleon-D-Ac coupling constant. 

We employ pointlike couplings between the current and the intermediate meson and 
baryon. For the vector current one has 

{Uv')\ciAUp)) = Wh.uip) (4) 

and 

(Dip')\c^^c\Dip)) = -{p + p'), (5) 

in a convention where the c-quark has charm charge=-|-l. 

The effective lagrangian eq.(|]) is non-local and this induces an electromagnetic vertex 
current if the photon is present. In order to maintain gauge invariance we have to take into 
account the "seagull vertex" 

Fik"^) - F((a ± kV) 
^r,{k, q) = ±g^^^T,^,{q ± 2k), , (6) 
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which is generated via minimal substitution . The upper and lower signs in eq. (6) corre 



spond to an incoming or outgoing meson respectively. 

The three distinct contributions to the intrinsic form factors, associated with processes 
in which the current couples to the baryon line (B) (figure la) , to the meson line (M) (figure 
lb) or to the meson-baryon vertex (V) (figure Ic and Id) in the loop are given by 



(2 

-A((k 

(27r 



X 



In the above equations 



(9) 



is the meson propagator and 



S{p,k)= \ (11) 

p — If — M\ + le 

is the Ac propagator and p' = p + q with q being the photon momentum. In figure 1 we show 
all momentum definitions. 

With these amplitudes it is easy to show that the Ward-Takahashi identity 

Q'Kip'.p) + r;'(p',p) + rj'(p',p)) = Qe(s(p) - s(pO) , (12) 

is satisfied. In eq.(12) Qc is the nucleon charm charge, Qc = 0, and S(p) is the self-energy 
of the nucleon related to the D loop. The sum of the three amplitudes also ensures the 
charge non-renormalization (or the Ward Identity) 

K + rf + r;:),=o = Qc (-^.m) = . (13) 



The intrinsic charm form factors are obtained by writing these amphtudes in terms of 
the Dirac and Pauh form factors 

T,{p',p) = l,F,%q') + ^^F^{q') . (14) 

The intrinsic squared charm radius of the nucleon is defined as 
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where C^^lq"^) is the electric form factor introduced by Sachs |TT 



(15) 



G'kiq') = F^{q') + ^F^{q') . (16) 



The numerical results for |r^| are shown in figure 2 , as a function of the form factor 
cut-off A. The value of the coupling and masses used are Mjy = 939 MeV, Ma^ = 2285 MeV, 
mjy = 1865 MeV and g^Ajj/V^ = Qn-kn/V^ = —3.795 

The intensity of a given proton fluctuation is associated with its average squared radius 
|r^|. The larger is \r^\ , the more frequently we will find the proton in that particular 
oscilation and the larger will be the probability of "seeing " it. 

We shall assume that the average barionic radius of the proton {vp = [< r\ >] ^ , ~ 
0.72 fm) associated with the isoscalar part of the electromagnetic current is a good measure 
of the proton "total size", i.e., the size which takes into account all possible fluctuations that 
couple to isoscalar currents. The intrinsic charm probability is then given by 



I T''^ 



P^c = ^ = 0.9% (13) 



where |r^| = 0.0047/m^ is the average charm squared radius calculated above with a cut-off 
A = 1.2 GeV. Pic is the ratio between the charm "area" and the total proton "area". 

We want to compare our results with those obtained by Donoghue and Golowich in ref. [5] 
for the five quark components of the proton wave function , \uudss > and \uudqq > , where 



q represents a light quark. We repeat then the calculations for kaon and pion loops ( with 
the same cut-off A) , obtaining the average strange radius \r'^\ = 0.025 fm'^ and the average 
light quark radius |r^| = 0.130/m^. Dividing these radii by the barionic squared radius used 
above we obtain the probabilities Pis = 5 % and Pig = 25 % . The calculations done in 
ref. [5] arrive at Pis = 16 % and Piq = 31 % . The discrepancy in the strange sector suggests 
that the vector meson dominance model contribution coming from the uo — cf) mixing ( see 
ref. [9]) is really important . In fact, it will change the result from Pis = 5 % to Pis = 10 
% [9]. As there is no experimental evidence for a — J/\I' mixing , the vector meson model 
will not contribute in the charm sector. With the inclusion of the u — (p mixing our results 
agree with those obtained in ref. [5] within 6 %. 

The charm squared radius increases with A (as it can be seen in figure 2) reaching 
|r^| = 0.016/m^ at asymptotically large values of A. In this limit we would have Pic = 3.0 
%. Considering that we are overestimating the coupling constant in the charm loop, this 
number can be taken as an upper limit for the intrinsic charm probability in the context of 
our calculation scheme. Our result seems to corroborate the previous estimates [1-5] . 

In this work we have only considered loops involving the particular combination D-A^ . 
In principle we could include loops with -D-Sc and also with vector mesons. However in the 
case of the intrinsic strangeness these contributions were shown to be much less important 
than the D-Ac loops. The couples very weakly to the nucleons and the vector mesons 
have large masses , their contribution being thus supressed. In our case, due to the lack of 
knowledge of the relevant couplings and cut-offs , no attempt is made to go beyond the D-Ac 
loop. We expect this contribution to be the most significant, specially in view of the very 
large values of the coupling constant and cut-off used here. This might be sufficient for an 
estimate of the order of magnitude of P^c- 
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Figure Captions 



Figure 1. Diagrams wliicli contribute to the calculation of the vertex function. Solid 
external lines represent the proton and solid internal lines represent the Ac. Dashed and 
wavy lines represent the D and the vector current respectively. 

Figure 2. The intrinsic charm mean square radius of the nucleon as a function of the cut-off 
A in the baryon-meson form factor. 
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